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We revisit the path integral description of the motion of 
a relativistic electron. Applying a minor but well motivated 
conceptional change to Feynman's chessboard model, we ob- 
tain exact solutions of the Dirac equation. The calculation is 
performed by means of a particular simple method different 
from both the combinatorial approach envisaged by Feynman 
and its Ising model correspondence. 



Here, N is the number of segments of constant length 
e = t/N of the particle's path between its start point 
(which is assumed to be the origin of the corresponding 
coordinate system) and the end point {x, t) of the path. 
R denotes the number of bends while Ns^ (R) stands for 
the total number of paths consisting of N segments with 
R bends. The indices 7 and S correspond to the directions 
forward or backward at the path's start and end points, 
respectively, and refer to the components of K. Asj{e) 
accounts for a convenient normalization. 
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I. INTRODUCTION 

It is well known that the continuum propagator of 
the Dirac equation can be found by summing over ran- 
dom walks. Renewed interest in this issue has arisen in 
connection with the investigation of stochastic processes 
which have been shown to be related to the Dirac equa- 
tion (Gaveau et al. 1984, McKeon and Ord 1992). Like- 
wise, the correspondence between the path integral and 
the Ising model has been explored (Gersch 1981, Jacob- 
son and Schulman 1984) and solutions for a discretized 
version of the Dirac equation have been found (Kauffman 
and Noyes, 1996). 

As described by Feynman and Hibbs (1965), the prop- 
agator of the 1 + 1 dimensional Dirac equation 



(1) 



(where units c = h = I are assumed and cTx and CTz are 
the respective Pauli spin matrices) can be found from a 
model of the one-dimcnsional motion of a relativistic par- 
ticle. In this model the motion of the electron is restricted 
to movements either forward or backward occuring at the 
speed of light. Assuming units c = 1, the motion of the 
particle corresponds to a sequence of straight path seg- 
ments of slope ±45° in the x-t plane. The retarded prop- 
agator K(x,t) of the Dirac equation may then obtained 
from the limiting process (see e.g. Feynman and Hibbs 
1965, Jacobson and Schulman 1984) 
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II. MODEL AND CALCULATIONS 

In this short note we demonstrate that a minor concep- 
tional change of Feynman's chessboard model naturally 
and directly yields exact solutions to the Dirac equation 
(|l|). The conceptional change is suggested by the obser- 
vation that a path with R bends between given start and 
end points is determined by _R — 1 bends. For a sketch of 
the situation consider Figure |^. The path shown in Fig- 
ure |l| exhibits five bends, three to the left and two to the 
right. However, the first two bends to the right and left, 
respectively, determine the path since the location and 
direction of the last bend (indicated by a circle in Figure 
|l|) is fully determined by the first four bends. We thus 
consider here, in contrast to the original formulation of 
the model where all bends occuring on a path contribute 
to the total amplitude, only contributions to the total 
amplitude from bends which actually define the path. In 
the light of the general path integral formalism it makes 
perfect sense to consider only those bends which define 
a path, i.e. the minimum information characterizing a 
path. 

In the following we demonstrate by an explicit calcu- 
lation that the modified model directly leads to exact 
solutions of the Dirac equation (|l|). We will use a calcu- 
lation scheme different from the combinatorial approach 
envisaged by Feynman and Hibbs (1965) and its Ising 
model correspondence (Gersch, 1981). Following Feyn- 
man's chessboard model we consider each bend which 
defines a possible path to contribute an amplitude 

(t>jr = (3) 

where ej^ = e is the length of a path segment. The total 
amplitude contributed by a path is the product 
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where jr runs over all the segments followed by a bend. 
While the index r enumerates the path segments after 
which bends occur, the value of jr indicates the corre- 
sponding segment. A path with R bends which starts 
with positive velocity (i.e. to the right) and ends with 
negative velocity (i.e. to the left) consists of exactly 
(i? - l)/2 + 1 bends to the left and {R - l)/2 to the 
right. The {R— l)/2 bends to the right may occur after 
any arbitrary path segment to the left, (i? — l)/2 of the 
(i? — l)/2 + 1 bends to the left occur in the same man- 
ner after path segments to the right while the additional 
bend to the left occurs after the last segment. Let P be 
the total number of path segments to the right (-1-) and 
Q those to the left (-). Then, the contribution of the 
i?+ = (i? - l)/2 bends to the right to ^'_+ is 



^^+{R+)=N_+{R+)Y[imie,^) 



jl<---<jR + 

For P 1, \1/ ^(i?^) is approximated by 

*_+(i?+) 
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The contribution of the i?" = [((i? - l)/2 + 1) - 1] = 
{R — l)/2 bends to the left is calculated similarly. The 
additional bend (occuring after the last segment to the 
right) does not enter the calculation since a possible path 
is fully determined by the location by its i? — 1 bends to 
the right and left, respectively. Therefore we find 
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In the limit N ^ oo (i.e. P,Q oo) the exact expression 
for 'J L becomes 
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where e = t/{P + Q). With v = Ax/At = x/t 
{P ~ Q)/{P + Q) the classical velocity attributed to the 
particle, PQ = ((P + Q)/2'yf where 7 
Thus we have 
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where Jq is the zeroth order Bessel function of the first 
kind. A similar calculation yields for the same re- 
sult. 

For 4"++, the number of bends to the right and to the 
left is i?/2 for each direction where R is even. However, 
the path is again defined by i?+ — R/2 bends to the right 
and R- = R/2 - 1 bends to the left. Thus, 
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(i?/2)!(P/2-l)! 
i^PlQ.h{mth) . 



(10) 



With ^/pJQ^ [t + x)l(f -x^ fl'^ andT= (t2-a;2)i/2 
the component 4*++ becomes 



*++ = i(t + a;)/r Ji(mt/7) . 
A similar calculation yields 

* = i(t - a;)/T Ji (mt/7) . 
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This completes the envisaged computation. As a side 
remark note that the presented calculation scheme is not 
restricted to tj^ = e. As will be shown elsewhere, similar 
results may be obtained for cj^ = e{jr)- 



III. DISCUSSION 

To relate the components 'i^s-y to the solution of the 
Dirac equation (^ consider the explicit represation 
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As may be seen by direct calculation, in this representa- 
tion "ifi and 5" 2 defined as 



*i = 



^9 = 



(14) 



are two independent, exact solutions of the Dirac equa- 
tion (|l|). This completes the demonstration that Feyn- 
man's chessboard model yields exact solutions to the 
Dirac equation when taking into account only those 
bends which actually define paths. With regard to funda- 
mental theories of spacetime and/or quantum mechanics 
(e.g. in the spirit of Finkelstein, 1974) this could be of 
importance. Somehow similar results have been obtained 
from the continuum limit of a discretized version of the 
Dirac equation (Kauffman and Noyes, 1996). 
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The calculation scheme and part of the results pre- 
sented here can be generalized to unevenly spaced space- 
time lattices. This opens up the possibility to de- 
fine an analogon to the Feynman checkerboard for dis- 
crete spacetime models of different type (e.g. KuU and 
Treumann, 1994). Related work is in progress and will 
be presented elsewhere. 
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FIG. 1. A possible path with five bends between given start 
and end points. The first two bends to the right and left, re- 
spectively, determine the path since the location and direction 
of the last bend (indicated by a circle) is fully determined by 
the first four bends. 
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